UNIVERSITY  OF  MARYLAND 
COMPUTER  SCIENCE  CENTER 


COLLEGE  PARK,  MARYLAND 


V - — 


"  ~  ^  3^  ^'Z>AKPA  <Sfd*r-~2JpC 


1 


TR-755  J 

DAAG-53-76C-0138 


Ap 


COMPUTING  _PER I METERS  OF  .IMAGES 
REPRESENTED  BY  QUADTREES* 

<  <Bir 

If)  I  Hanan/samet  . 
CompinSerScience  Department 
University  of  Maryland 
College  Park,  MD  20742 


(ffiirn?ch 


79  / 

-J 


n  D  Q, , 

\'J'  offi  '  1919  \\ 


ll;.  •  > 
U  UD:- 


ABSTRACT 

. 

An  algorithm  is  presented  for  computing  the  total 
perimeter  of  a  binary  image  represented  by  a  quadtree. 
The  algorithm  explores  each  segment  of  the  border  once 
and  only  once.  Analysis  of  the  algorithm  shows  that 
its  worst-case  average  execution  time  is  proportional 
to  the  product  of  the  log  of  the  image  diameter  and 
the  number  of  nodes  in  the  tree. 
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1.  Introduction 


Perimeter  computation  i  ;  a  basic  operation  in  image  proces¬ 
sing  IRK] .  The  standard  algorithms  use  either  an  array  or  a 
chain  code  representation  [:  reeman]  for  the  two-valued  ("binary") 
image  whose  perimeter  is  to  be  computed.  In  this  paper  we  present 
an  algorithm  for  computing  the  total  perimeter  (i.e.,  the  length 
of  the  chain  codes  corresponding  to  the  black/white  borders  in  the 
image)  of  a  binary  image  that  is  represented  by  a  quadtree 
([Klinger,  DRS ,  Sametl] ) . 

We  assume  that  the  given  image  is  a  2n  by  2n  array  of  unit 
square  "pixels".  The  quadtree  is  an  approach  to  image  repre¬ 
sentation  based  on  successive  subdivision  of  the  array  into 
quadrants.  In  essence,  we  repeatedly  subdivide  the  array  into 
quadrants,  subquadrants,...  until  we  obtain  blocks  (possibly 
single  pixels)  which  consist  entirely  of  either  l's  or  0's.  This 
process  is  represented  by  a  tree  of  out-degree  4  in  which  the 
root  node  represents  the  entire  array,  the  four  sons  of  the  root 
node  represent  the  quadrants,  and  the  terminal  nodes  correspond 
to  those  blocks  of  the  array  for  which  no  further  subdivision  is 
necessary.  For  example.  Figure  lb  is  a  block  decomposition  of  the 
region  in  Figure  la  while  Figure  lc  is  the  corresponding  quad¬ 
tree.  In  general,  E’LACK  anc  WHITE  square  nodes  represent  r.odes 
consisting  entirely  of  l's  and  0's  respectively.  Circular  nodes, 
also  termed  GRAY  nodes,  denote  non-terminal  nodes. 


Sections  2-5  present  and  analyze  the  algorithm.  Included 
is  an  informal  description  of  the  algorithm  along  with  moti¬ 
vating  considerations.  The  actual  algorithm  is  given  using  a 
variant  of  ALGOL  60  [Naur] . 
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2 .  Definitions  and  Notation 

Let  each  node  in  a  quadtree  be  stored  as  a  record  contain¬ 
ing  six  fields.  The  first  five  fields  contain  pointers  to  the 
node's  father  and  its  four  sons,  labeled  NW,  NE,  SE,  and  SW. 
Given  a  node  P  and  a  son  I,  these  fields  are  referenced  as 
FATHER (P.)  and  SON(P,I)  respectively.  At  times  it  is  useful  to 
use  the  function  SONTYPE(P)  where  SONTYPE(P)  =  Q  iff  SON (FATHER 
(P) ,Q)  =  P.  The  sixth  field,  named  NODETYPE,  describes  the 
contents  of  the  block  of  the  image  which  the  node  represents — 
i.e.,  WHITE  if  the  block  contains  no  l's,  BLACK  if  the  block 
contains  only  l's,  and  GRAY  if  it  contains  0's  and  l's.  Alter¬ 
natively,  BLACK  and  WHITE  are  terminal  nodes  while  GRAY  nodes 
are  non-terminal  nodes . 

Let  the  four  sides  of  a  node's  block  be  called  its  N,  E,  S, 
and  W  sides.  They  are  also  termed  its  boundaries.  The  spatial 
relationships  between  the  various  sides  are  specified  by  use  of 
the  functions  OPSIDE,  CSIDE,  and  CCSIDE.  OPSIDE(B)  is  a  side 
facing  side  B;  e.g.,  OPSIDE (E)  =  W.  CSIDE (B)  and  CCSIDE (B)  cor¬ 
respond  to  the  sides  adjacent  to  side  B  in  the  clockwise  and 
counterclockwise  directions  respectively;  e.g.,  CSIDE(E)  =  S  and 
CCSIDE (E)  =  N.  We  also  define  the  following  predicates  and 
functions  to  facilitate  the  expression  of  operations  involving  a 
block's  quadrants  and  boundaries.  ADJ(B,I)  is  true  if  and  only 


if  quadrant  I  is  adjacent  to  boundary  B  of  the  node's  block; 
e.g. ,  AD J (N, NW )  is  true.  REFLECT (B, I)  yields  the  quadrant 
which  is  adjacent  to  quadrant  I  along  boundary  B  of  the  block 
represented  by  I;  e.g.,  REFLECT (W,NW)  =  NE,  REFLECT (E ,NW)  =  NE, 
REFLECT (N , NW)  =  SW,  and  REFLECT ( S ,NW)  =  SW.  QUAD(B,C)  is  the 
quadrant  which  is  bounded  by  boundaries  B  and  C  (if  B  and  C  are 
not  adjacent  boundaries,  then  the  value  of  QUAD(B,C)  is  undefined) 
e.g.,  QUAD (N , W)  =  NW.  Figure  2  shows  the  relationship  between  the 
quadrants  of  a  node  and  its  boundaries. 

Given  a  quadtree  corresponding  to  a  2n  by  2n  array,  we  say 
that  the  root  node  is  at  level  n,  and  that  a  node  at  level  i  is 
at  a  distance  of  n-i  from  the  root  of  the  tree.  In  other  words, 
for  a  node  at  level  i,  we  must  ascend  n-i  FATHER  links  to  reach 
the  root  of  the  tree.  Note  that  the  farthest  node  from  the  root 
of  the  tree  is  at  level  £0.  A  node  at  level  0  corresponds  to  a 
single  pixel  in  the  image. 
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Figure  2.  Relationship  between  a  block's  four  quadrants 
and  its  boundaries. 


3 .  Algorithm 


The  perimeter  computation  algorithm  traverses  the  quadtree 
in  postorder  (i.e.,  the  sons  of  a  node  are  visited  first)  and 
visits  each  BLACK-WHITE  border  segment  once  and  only  once. 

For  each  BLACK  terminal  node,  say  P,  the  northern,  eastern, 
southern,  and  western  adjacencies  are  explored  so  that  all  of 
the  node's  WHITE  adjacent  neighbors  are  visited.  The  result  of 
each  visit  is  that  the  length  of  the  border  which  is  shared 
between  the  two  adjacent  neighbors  is  included  in  the  value  of 
the  perimeter. 

The  main  procedure  is  termed  PERIMETER  and  is  invoked  with 
a  pointer  to  the  root  of  the  quadtree  representing  the  image 
and  an  integer  corresponding  to  the  log  of  the  diameter  of  the 
image  (e.g.,  n  for  a  2n  by  2n  image  array).  PERIMETER  traverses 
the  tree  and  controls  the  exploration  of  the  adjacencies  of  each 
BLACK  node.  FIND_NEIGHBOR  locates  a  neighboring  node  of  greater 
or  equal  size  along  a  specified  side.  If  the  node  is  on  the 
edge  of  the  image,  then  no  neighbor  exists  in  the  specified 
direction  and  NULL  is  returned  (e.g.,  border  segments  CO  and  CN 
in  Figure  lb) .  If  the  node  is  not  on  the  edge  of  the  image  and  nr 
neighboring  BLACK  or  WHITE  node  exists  satisfying  our  size  cri¬ 
teria,  then  a  pointer  to  a  GRAY  node  of  equal  size  is  returned 
(e.g.,  the  eastern  border  of  node  C  in  Figure  lb).  In  such  a  c<.»>:, 
procedure  SUM_ADJACENT  continues  the  search  by  examining  all  WHITE 


adjacent,  neighbors  of  smaller  size  and  accumulating  their  lenqths 
(e.g.,  block  M  for  the  eastern  border  of  node  C  in  Figure  lb). 
Note  that  when  the  neighboring  node  is  BLACK  and  is  of  the  same 
or  greater  size,  then  no  contribution  is  made  to  the  perimeter 
by  the  side  of  the  node  currently  being  examined  (e.g.,  border 
segment  AB  in  Figure  lb) .  A  node  having  a  side  on  the  border 
of  the  image  or  having  a  WHITE  neighboring  node  of  the  same  or 
greater  size  makes  a  contribution  to  the  perimeter  equal  to  the 
length  of  the  side  of  the  BLACK  node  (e.g.,  border  segment  AI 
in  Figure  lb) . 

An  alternative  method  of  computing  the  perimeter  is  to  apply 
the  algorithm  in  [DRS]  which  converts  a  quadtree  representation 
to  a  chain  code  and  simply  sums  the  lengths  of  the  segments.  Our 
algorithm  is  simpler  since  it  does  not  require  the  segments  to 
be  traversed  in  sequence  around  each  border.  We  need  only  in¬ 
sure  that  each  border  segment  is  visited  once  and  only  once. 

This  is  clearly  true  since  during  the  tree  traversal,  the  adja¬ 
cencies  of  each  BLACK  node  are  explored  at  least  once;  on  the 
other  hand,  each  border  segment  is  only  explored  once  since  it 
must  adjoin  a  WHITE  node  and  our  algorithm  does  not  explore  adja¬ 
cencies  of  WHITE  nodes. 

As  an  example  of  the  application  of  the  algorithm,  consider 
the  image  given  in  Figure  la.  Figure  lb  is  the  corresponding 
block  decomposition  and  Figure  lc  is  its  quadtree  representation 
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All  of  t  ie  BLACK  nodes  have  labels  rangi  ng  between  A  and  G 
while  the  WHITE  nodes  have  labels  ranging  between  H  and  S.  The 
BLACK  nodes  are  labeled  in  the  order  in  which  their  adjacencies 
are  explored  by  PERIMETER.  WHITE  nodes  H  through  Q  are  labeled 
in  the  order  in  which  they  are  first  visited  by  the  combination 
of  FIND_NEIGHBOR  and  SUM_ADJACENT .  Thus  the  adjacencies  of 
node  A  have  been  explored  before  those  of  nodes  B,  C,  etc.  The 
value  of  the  perimeter  is  obtained  by  v. siting  the  border  seg¬ 
ments  in  the  order  AH,  AI ,  BJ,  BK,  BL,  Cl,  CM,  CN,  CO,  EL,  EP , 
FP,  GQ,  and  GM.  Assuming  n=3  (i.e.,  blocks  D,  E,  F,  G,  P,  Q,  R, 
and  M  are  single  pixels),  the  perimeter  is  28.  Note  that  nodes 
D,  R,  and  S  do  not  contribute  to  the  value  of  the  perimeter 
since  none  of  their  sides  adjoin  the  border. 


integer  procedure  PE RIMETER(P, LEVEL) ; 

/*find  the  perimeter  of  a  quadtree  rooted  at  node  P  which 
spans  a  2LEVEL  by  2LEVEL  space*/ 
begin 

node  P,Q; 

integer  LEN ,  LEVEL ; 
quadrant  I; 
side  S; 

LEN<-  0 ; 

if  GRAY (P)  then 
begin 

for  I  in  {  NW,NE,SW,SE)  do 

LEN+-LEN+PERIMETER  ( SON  ( P  ,  I )  ,  LEVEL- 1)  ; 


end 

else  if  BLACK (P)  then 
begin 

for  S  in  {N,E,S,wl  do 


begin 


Q<-FIND_NEIGHBOR(P,S)  ; 

LENGLEN  +  if  NULL (Q)  o_r  WHITE  (Q)  then  2+  LEVEL 
else  if  GRAY (Q)  then 

SUM_ADJACENT (Q,QUAD (OP SIDE (S) ) ,CSIDF(S) 

OUAD (OPSIDE (3) ,CCSIDE(S) ) 
LEVEj  ) 

else  0; 


end; 

end; 

return (LEN) ; 


k 


end; 


in  direction  S.  Of  these  2n  1-1  neighbor  pairs,  2^  have  their 
nearest  common  ancestor  at  level  n,  2*  at  level  n-1,...  and 
2n  1  1  at  level  i+1.  For  each  node  at  level  i  having  a  common 
ancestor  at  level  j,  the  maximum  number  of  nodes  that  will  be 
visited  by  FIND_NEIGHBOR  and  SUM_ADJACENT  is 

(j-i)  +  ( j-i-1)  +  £  2k  =  2 ( j-i-1)  +  21+1 

k=0 

This  is  obtained  by  observing  that  the  common  ancestor  is  at 
a  distance  of  j-i  and  that  a  node  at  level  i  has  a  maximum  of 
21  adjacent  neighbors  (all  appearing  at  level  0) .  Assuming  that 
node  P  is  equally  likely  to  occur  at  any  level  i  and  at  any  of 
the  2n  1-1  positions  at  level  i,  then  the  average  of  the  maximum 
number  of  nodes  visited  by  FIND_NEIGHBOR  and  SUM_ADJACENT  is 
n-l  n  . 

£  £  2  3 (2 (j-i-1) +2  ) 

i=0  j  =  i+l _  (1) 

n  _  j 

£  ( 2n  -1) 
i=0 

(1)  can  be  rewritten  to  yield 

n-l  n-l-i  ,  .  .  . ,, 

£  £  2  1  1-3 (2 j+21+1) 

i=0  j=0 _ 

n  i 

£  (2-1) 

i=0 


(2) 


The  numerator  of  (2)  can  be  simplified  as  follows 


n-1  n-l-i  ,  .  . 

I  E  2n  ^  1  ^ 
i=0  j=0 


(2j+21+1) 


n-1  n-l-i 

E 

i=0  j=0 


E  (j-2n_1~j+2n 


n-1  n-l-i  /  . 

=  E  2n  Z 

i=0  j=0  '211 


n-l-i 


i+D 


But  n-l-i 

E 


,  n-l-i 

i .  4  i  4 

j_0  21+-^  21  j=o  2^ 


-i.(2 _ n-i-, ) 

2n-2-i) 


Also  n-l-i  -i 

E  A  =  2(1  “  ^ 

j=0  2J  2 


Substituting  (4)  and  (5)  into  (3)  yields 

n"1!n  (  -4(2  -  ^T>  +  2U  -  2. 


n-1  „  /  ■■ 

£  2  (4 

i=0  \2 


-pM 


n-1  /-n-l-i  ,  ..  -n+l-i  „\ 

=  Z2n(? — .  4.  2 .  -f--2) 

i=0  1  2n'2  2  1  ' 


n-1 


=  E  (2n+1-1-4 (n-i)+2n+1-21+1) 


i=0 


,  i  n~l  i  n-1  1 1  n  1  •  .  -1 

=  2n  1  E  — i-  -  4  E  (n-i)  +  n-2  1  -  E  2  1 
i=0  21  i=0  i=0 


=  2°+1  (2(1 - — )  )  -  4n--j}+ +  n •  2n+1-2n+1+2 

2  n  4 


=  2n+2  -  4  -  2n2  -  2n  +  (n-1) 2n+1  +  2 


=  2n+1 (n+1)  -  2 (n2+n+l) 


4 .  Analysis 

The  running  time  of  the  perimeter  computation  algorithm, 
measured  by  the  number  of  nodes  visited,  depends  on  the  time 
spent  locating  adjacent  WHITE  nodes  and  on  the  size  of  the  quad¬ 
tree.  Adjacent  WHITE  nodes  are  located  by  procedures  FIND_NEIGHBOR 
and  SUM_ADJACENT .  They  are  invoked  once  for  each  BLACK  node. 

The  amount  of  work  performed  by  these  procedures  is  obtained  by 
considering  the  number  of  nodes  that  are  visited  when  an  adja¬ 
cency  is  being  explored.  Recall  that  we  must  find  the  neighbor, 
and  if  it  is  GRAY,  then  visit  all  adjacent  WHITE  neighbors  of 
smaller  size.  In  the  worst  case  we  are  at  level  n-1,  with  a 
GRAY  neighbor,  and  all  adjacent  neighbors  are  at  level  0.  In  such 
a  case,  we  must  visit  2n  nodes.  For  example,  consider  Figure  lb 
where  n=3  and  we  wish  to  visit  the  blocks  adjacent  to  the  block 
labeled  C  (i.e. ,  blocks  D,  F,  G,  and  M) .  We  must  visit  the  root 
of  the  quadtree  as  well  as  A's  neighboring  GRAY  node  and  a .1 1  of 
its  NW  and  SW  sons — i.e.,  a  complete  binary  tree  of  height  2. 

In  total,  2^=8  nodes  are  visited.  Assuming  a  2n  by  2n  random 
image--i.e.,  a  BLACK  node  is  equally  likely  to  appear  in  any 
position  and  level  in  a  quadtree,  we  have  the  following  result. 

Theorem  1 :  The  average  of  the  maximum  number  of  nodes  vis.’ ted 
by  each  invocation  of  FIND  NEIGHBOR  and  SUM_ADJACENT  is  n+] . 

Proof :  Given  a  node  P  at  level  .i  and  a  direction  S,  there  art 
2n-1-l  possible  positions  for  node  P  and  a  neighbor  at  level  i 


node  procedure  FIND_NEIGHBOR(P,S) ; 

/♦given  node  P,  return  a  node  which  is  adjacent  to  side  S  of  node  P*/ 
begin 

node  P,Q; 
side  S; 

if  not  NULL (FATHER (P) )  and  ADJ (S,SONTYPE (P) )  then 
/*  find  a  common  ancestor*  / 

Q<-FIND_NEIGHBOR (FATHER (P)  ,S) 
else  Q-*-FATHER (P)  ; 

/♦follow  reflected  path  to  locate  the  neighbor*/ 

return  (if  not  NULL (Q)  and  GRAY (Q)  then  SON (Q, REFLECT (S ,SONTYPE P) ) ) 
else  Q) ; 

end; 

integer  procedure  SUM_ADJACENT (P ,Q1 ,Q2 , LEVEL) ; 

/♦find  all  WHITE  descendants  of  node  P  adjacent  to  the  perimeter - 
i.e.,  in  quadrants  Q1  and  Q2,  and  return  the  length  of  their  ~ 
begin 

node  P; 

quadrant  Q1,Q2 
integer  LEVEL; 

return  (if  GRAY(P)  then  SUM  ADJACENT (SON (P ,Q1)  . Q1 ,Q2 ,LEv £L- 1 i 

SUM_AD JACENT  ( SON  ( P , Q2 )  , Q1 , Q2  , LEV L..-1 
else  if  WHITE  (P)  then  2+LEVEL 
else  0) ; 


end; 


The  denominator  of  (2)  can  be  simplified  as  follows: 


?  ( 2 1— 1 )  =  2n+1-l- (n+1) 
i=0 


or  n 


E  (21-!)  =  2n+1-n-2 


i=0 


Substituting  (6)  and  (7)  into  (2)  yields 

2n+1 (n+1) -2 (n2+n+l)  =  n+1  _  n2-n 

-n+1  -  -n+1  0 

2  -n-2  2  -n-2 


(7) 


=  n+1  as  n  gets  large 
Q .  E .  D . 

It  is  useful  to  obtain  an  upper  bound  on  the  size  of  the  quad¬ 
tree  in  terms  of  the  number  of  BLACK  nodes.  Letting  B  denote 
the  number  of  BLACK  nodes  we  have 


Lemma  1 :  The  maximum  number  of  nodes  in  a  quadtree  having  B 
black  nodes  is  4Bn+l. 

Proof:  Given  B  BLACK  nodes  at  level  0  there  is  a  maximum  of 

3B  WHITE  nodes  at  level  0.  At  worst  there  is  one  GRAY  node  at 
level  1  for  every  BLACK  node  and  three  WHITE  nodes  at  level  0, 
and  at  worst  three  additional  WHITE  nodes  for  each  such  GRAY 
node.  Thus  there  exist  at  most  4B  nodes  at  level  1.  Repeating 
the  same  argument  for  levels  2  through  n-1,  we  have  4Bn  nodes. 
At  level  n  there  is  only  one  node.  Therefore,  the  maximum 
number  of  nodes  is  4Bn+l. 


Q.E.D. 


We  can  now  prove 

Theorem  2:  The  average  worst  case  execution  time  of  the 
perimeter  computation  algorithm  has  an  upper  bound  proportional 
to  the  product  of  the  number  of  BLACK  nodes  and  the  log  of  the 
diameter  of  the  image. 

Proof :  From  Theorem  1  we  have  that  for  each  adjacency  involving 
a  BLACK  node,  FIND_NEIGHBOR  and  SUM_ADJACENT  result  in  an 
average  worst  case  of  n+1  nodes  being  visited.  There  are  four 
adjacencies  for  each  BLACK  node.  Thus  these  two  procedures  con¬ 
tribute  4B(n+l).  From  Lemma  1  we  have  that  the  number  of  nodes 
in  the  quadtree  is  bounded  by  4Bn+l.  However,  this  quantity  cor¬ 
relates  with  the  work  performed  by  procedure  PERIMETER  since 
each  node  in  the  quadtree  is  visited  by  the  traversal.  Summing 
up  these  values  we  have  4B(n+l)  +  4Bn+l  =  8Bn+4B+l. 


Q.E.D. 
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5.  Concluding  remarks 


An  algorithm  has  been  presented  for  computing  the  total 
perimeter  of  a  binary  image  represented  by  a  quadtree.  The 
algorithm's  execution  time  has  been  shown  to  have  an  average 
worst  case  time  complexity  proportional  to  the  product  of  the 
image's  diameter  and  the  number  of  BLACK  nodes  in  the  quadtree 
representation  of  the  image.  It  should  be  clear  that  if  the 
image  has  more  than  one  connected  component,  the  algorithm 
returns  the  total  perimeter  of  all  the  regions.  Similarly,  if 
holes  are  present,  their  boundaries  are  also  included  in  the 
value  of  the  perimeter  obtained  by  this  algorithm.  Note  that 
if  we  first  labeled  the  connected  components  of  the  image 
[Samet2] ,  then  the  perimeter  o f  each  boundary  could  be  separ¬ 
ately  computed. 

The  algorithm  demonstrates  the  utility  of  the  quadtree  as 
a  desirable  data  structure  for  image  representation.  Computation 
of  perimeter  is  generally  achieved  by  use  of  a  chain  code  repre¬ 
sentation.  We  have  shown  that  it  can  be  computed  with  reasonable 
efficiency  when  the  quadtree  is  used  as  the  data  structure. 
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An  algorithm  is  presented  for  computing  the  total  perimeter  of  a 
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».  SAMPLE  IMAGE  b.  BLOCK  DECOMPOSITION  OF  THE  IMAGE 
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c.  QUADTREE  REPRESENTATION  OF  THE  BLOCKS 
IN  (b). 

FIG.  1.  AN  IMAGE,  ITS  MAXIMAL  BLOCKS,  AND  THE  CORRESPONDING  QUADTREE.  BLOCKS 
IN  THE  IMAGE  ARE  SHADED. 


